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O . Abstract 
(N 

^ ' After the gauge invariant gluon condensation, gluons remain as massless excitations. When an effective 

theory desribing the condensation and the excitaion is constructed, a constraint has to be imposed for the 
, vacuum to be stable. The constraint impHes the perfect dia-electricity and assures the sohition of color 

flux tube when qurks are introduced. Thus the dia-electric model of Kogut-Susskind 5 and 't Hooft [6j 
are derived by the stability of the condensed vacuum. 



Oh: 



PACS nimber; 11.15-q, 12.38Aw, 12.38Lg. 



1^ . Introduction; 

\G , Understanding of the confinement mechanism based on QCD Lagrangian remains as a fundamental 

00 ■ 

' problem in the strong interaction regime. Besides lattice QCD jlj, active researches have been performed 

' along the line of the dual superconductivity with the abelian gauge fixing [51 [3] , infrared structure of the 

o ; 

00 ! gluon propagators [4], e.t.c. In continuum QCD, the key point will be how to extract the effect of the 

gluon condensation. 

, , In this letter, we propose an approach to construct the effective Lagrangian and derive the dielectric 

^ ' model of Kogut-Susskind [5] and 't Hooft [6J. We asssume that gluons condense in gauge invariant 

^ ' 

form, so the gauge invariance is not broken in the vacuum and gluons which are excitations above the 
condensed vacuum remain massless. Because of masslessness, they may condense further in color singlet 
form. Therefore, when we construct the effective Lagrangian describing the interaction between the 
condensation and excitation field, non-trivial conditions have to be imposed to assure the the stability of 
the vaccum. Remarkably, these conditions imply the perfect dia-electricty of the vacuum and guarantee 
the existence of the color flux tube of infinite length when quarks are introduced. As for the gluon 
condensation, the magneteic type condensation <G^j^>> has been established [7j but the advantage of 
the present work is that we do not have to ask the precise mechanism of the condensation, the magnetic 
monopolejS] or the gluon pairs|9j e.t.c. 

Effective Lagrangian Ces; 

To measure the singlet condensation, any gauge invariant operator can be used but we choose for 
simplicity the normal oredred form of Gj^^; (j){x) = A/'G^^,(x) = Gfj^^{x)~ <0|G'^^(x)|0>. (|0> is the 



normal vacuum and the dimensional regularization is adopted.) With |0>c the condensed vacuum, we 
know that c<0|A/'G'^y |0>c= 0c > 0. We regard gluon fields A'^{x) and as independent fields. (The 
general formalism to achieve this is shown below. ) As in the case of the free energy of the phenomenological 
theory of Ginzburg-Landau, the low energy effective Lagrangian CcS is obtained by the hydrodynamic 
expansion. It agrees with the expansion according to the operator dimension in By gauge 

invariance and keeping lowest non-trivial terms, we get 

= 5^0(a;)a^0(x)/2 - V{cl>{x)) 

-e{^{x))G%{x)/4^C^+C,,A, (1) 
G^,(a:) = a^A^(x) - d,A';^{x)+gr'^^Al{x)A:{x). 

Here, g, f"''"^ is the coupling constant and the structure constant of QCD, respectively. The effect of the 
condensation on the excitation appears as a dielectric factor e{(j)). Assuming ipc > 0, define V^(0) = 0, 
V{(t)c) < 0. We require (j) = {(p = (pc) is the maximum (minimum) of V{(f)); V'{0) = V'{(t)c) = 0, 
V"{0) < 0, V"{(t>c) > 0. When = 0, Cqcd has to be recovered, so e(0) = 1. For all <p, e(0) has to be 
positive in order for C^^a to be a sensible theory. The gauge fixing term and resulting ghost fields are not 
written excplicitly. Our conclusions will not change if these are introduced since only a gauge invariant 
degree (p{x) is added to the ordinary QCD, where we know how to extract physical sectors [TO] . 

Stability of the condensed vacuum; 

Below, the fiuctuation of the condensation is neglected, so (j){x) is treated as a c-number field. In order 
for £cff to describe a consistent theory, i.e. <j) — <j)c actually realizes the lowest energy state of £cff , we have 
to require first of all that excitation fields AJ^(a;) do not condense any more when ~ (pc- If e(</>c) 0, 
Ce^A is proportional to Cqcd, so A°(a;) condenses in the color singlet channel. In that case has to 
be redefined, and after that e(0c) = is satisfied. Such a perfect dia-electricity of the stable vacuum is 
specific to QCD; the gauge invariance is preserved after the condensation. Note that if the excitation 
acquires a mass gap, as in the case of the superconductor, such an instability never occurs. Assuming 
the homogeneous case — constant, we require that the energy is indeed increased for (f> ^ (f>c- Consider 
the trace of the energy momentum tensor in n-dimension by using Ccs; 6^ = 4V^(0) + e{(j)){n — 4)G^^. 
In calculating the anomalous term for n — > 4, we have to define A/'G^^,, by subtracting the expectation 
value taken by the perturbative vacuum of the ordinary QCD. In this conection, notice that our origin 
of the energy is that of the perturbative vacuum of QCD, for example we have defined V{4) = 0) = 0. 
Thus the situation is that the whole system evolves with while the normal order has to be defined 
by the perturbative vacuum of Cqcd- In such a case, it is convenient to adopt the interaction picture 
with C^=i^A — Cqcd as the "free Lagrangian", and the rest as the "interaction part", the latter being 
converted into the state vector. We can minimize the expectation value of 0^ in the interaction picture 
thus defined. (See below for details. ) Now the operator evolves by Cqcd, so the trace anomaly of QCD 
holds; 

e{:(a;) = W{cj>) + e{cj,mg)/2g)jVG^J{x), (2) 
Here P{g) is the usual (3 function with g the renormalized coupling constant. We know that the excitation 



field Af, condenses with tlie amount 

B = {f3{g)/8g),<0\UG-^^{x)\0>, < 0. (3) 

With fixed </>, the energy density (measured from the normal vacuum) of such a state is given by e = 
V{4i) + e{(j))B. Here we have used the relation valid for any homogeneous vacuum; <0^y(x)>= 5piye, so 
< >— Ae. We require the above condensation does not lower the energy of the state (f> — (pc, otherwise 
|0>c is unstable. Here we note that B introduced in ^ is the energy of condensed vacuum |0 >c, since 
it is calculated by Cqcd, therefore it is nothing but V{(t)c)', l^(0c) = B. In this way, the stability of |0 >c 
requires 

B < V{cl)) + e{4,)B, (4) 

A remarkable fact is that the stability condition ([4]) assures the existence of the solution of color flux tube 
of infinite length when quarks are introduced. Near ~ 0c, let us put (j) — (pc-'r A0 in (j4]) and expanding 
V{4)) and e((/)) up to (Ac/))^. One arrives at 

e(0c) = e'{cP,) = 0, V"{cj),) + e"{cP,)B > (5) 

In case e((/)) is parametrized as e((/)) — C{(j)~ (pc)'^"' with C > near (j) = (pc, then a > 1. Expanding in cj) 
for small (p, one obtains from ^ 

e'(0) = 0, y"(0) + e"(0)B>0. (6) 

By V"{0) < and B < 0, we get e"(0) < 0. Thus near = 0, e((/>) behaves as 1 + a^^ with a < 0. 

The inequality ^ becomes an equality for (p = (pc and for = 0. This should be the case since both 
represent |0>c; for = the excitation condenses producing B. However we have to select = 0c as 
|0>c, because our starting definition of is that it represents the exitation field without the condesned 
part. If = is selected, the role of and A/'G^^ is interchanged, so we have to restart our discussions. 
We conclude that = 0c is the lowest state of energy for < < 0c, without A/'G^^, condensing anymore, 
since the condensation does not lower the energy. When = 0c, we have no observable effect of JVGf^^ 
due to e = 0. (See below for more details.) Thus our effective Lagrangian is (H]), with V{(j)) and e(0) 
satisfying Before showing the tube-like solution, we discuss two subjects whose results have been 
utilized above. 

Constructing C^s', 

We present first the method of obtaining Ccs by the Fourier tranform and its inverse, with and A^^ 
being treated as independent fields. Let ^[AJJ] be an arbitrary functional of A'^{x) and start from 

Z = I [dA^MA^] expi I LgcDd^x (7) 

Let 0(A'J^{x)) be a gauge invariant operator and insert a functional identity J[d(j)] J[dJ] exp i / J(a;)(0(ylj^(a;))- 
ip{x))d*x and write ^'[^J^] by a functional Fourier transform Then the integration over A° is done 

with Cqcd + j^A"''^ + JO{A'^) in the exponential. Writing the result as expiiy[j'', J], 

^ = / [dJ] J [df^] J m * [j-] exp iW [J, cP, j] 

{dAi\\{dmAiwvn4>,K\- (8) 



is obtained. Here W[J,(l),j] — W[J,j^] — J J{x)(j){x)d'^x. In arriving at ([5]), ^'[j^] has been transformed 
back to 5" [A"^] , and we have done integrations over and J. Now after ip integration, we get back to Cqcd 
since only identical transformations are done. However, when the stationary phase of (j) integration is 
present, it corresponds to a nontrivial phase. Note that if the integral defining Z is dominated by a special 
value of (j), Z cannot be transformed back to Cqcd, inequivalent vacuum being selected. Since ^'[^^] is 
arbitrary, T[(p, A'^] can be regarded as the Lagrangian of the condensed phase. (^'[A^] can be generalized 
to ^'[^Ji, 0] in ([H]).) We have checked our useful formula ([5]) for several model Lagrangians. In the case of 
QCD, we know that (p condenses, so T[(j), A^j] has a a non-trivial stationary solution corresponding to this 
condensate. Eq.((T]) is the hydrodynamic expansion of r[(/), A"] and we have regarded as a c-number field 
since (j) loses fluctuations by the stationary phase mechanism. The stationary condition is represented by 
the stationary equation of and we require Ci^^a not to cause the instability of the of V'{(f>) — 0. 

Normal ordering and interaction picture; 

Consider the expectation value of an arbitrary local operator Q{A'^^{x)) = Q{x) in the theory £ofi' for 
the homogeneous case of 4>{x) — (p. Omitting the c-number C^, here we introduce the U operator in the 
path integral form by 



U{t2,tx : e) = 




([dA^] contains the factor e(0).) The matrix element of U between the states |A^j(ti_2, ^c) > (in 
coordinate representation) is equal to the right-hand side with AJ^(x) fixed to the value Aj^(ii^2,aj) at 
ti^2- Now we want to rewrite U{t2,ti : e) by the interction representation where "free part" is defined 
to be e = 1. (This method of introducing the interaction representation in path-integral form is quite 
convenient since it keeps the Lorentz invariance explicitly.) By the homogeneity of the ststem, we take 
a; = 0. With |5'> an arbitrary state and taking the coordinate representation in mind, let us consider 

<Q(0)>=<*|C/t(0,-c» : e)g(0)[/(0,-c» : e)|^'> 

=<^,(f)\U\-oo,0: l)Q(0)C/(0,-oo ; l)|«',(/)>, 
|^',(^>= [/(-oo,0 : l)f/(0,-(X) : e)|*> . 

In this representation, Q evolves by Cqcd so for Q{0) = {n — 4)6"°^ (0), the ordinary trace anomaly ^ 
holds. Here the normal order is defined with respect to the perturbative vacuum of QCD. As for the 
state vector, if we choose j^*, (p >— \0>c, then it is the lowest energy state in the theory Cqcd obtaining 
®. 

A comment here; if we insert the adiabatic factor connecting Cqcd and CcS, the difference of the 
energy of the ground states of two theories (differnece of both-hand sides of (jj])) is given by the adiabatic 
formula obtained in the interaction representation TlJ. 

Color flux tube; 

When color sources are present, they couple to the excitation field A'^{x). We show here that the 
c-number solution of the infinite tube exists, where abelian components Aq and Aq are non-zero. This 



K{cl>,e)=v{<j>)-e{cl>)8 




Figure 1: Schematic behavior of K{(j), 
K{0,e) = -9. 



< 0(0) < <j>c and K{(j){Q),0) < -1 for 9 > 1. Note that 



is sufficient for the confinement of color non-singlet state. Such a solution, if it exists, has a large c- 
number value, so it is stable against quantum fluctuations of other giuonic fields. (Quark pair creation 
is not discussed.) Below, the color index a takes 3 or 8. Now consider the case of static quark-antiquark 
point-like source with the strength ±g{X'^ /2) at z — ±oo. (A° is the representation matrix of the quark.) 
Define G^q = e{<j))E^ {k ~ x, y, z), then div{e{(j))E°') is the color density. For an infinite tube, cj) depends 
only on p in the cylindrical coordinate {z,p,(p) and E'^ is directed along z axis, which is constant over 
all space beause of rot£^" — 0. The equation of (p is 



(9) 



where we have written — E^X"- /2, E"^ = i?" — El/i for any color of quark. We have also defined 

K{<ly,9)^v{4,)-e{4,)9, 9 = E^ /2\B\. 

with !;((/)) = V{(I))/\B\. Note that Ki^cj), 1) > -1 by (g]). The behavior of K{(j),9) is shown in the Figure. 
The solution of © with (/)(p) approaching (pc sX p — co is given by 4>{p) ^ 4>c + Dieyij>{~\/W p) + 
D2ey.v{VWp), with W = \B\K"{(j)c). We see > by Q if a > 1. (When a = I, E < E^^^ 
with El^^^ = 2\B\v"{(l)c)/e"{(t>c) has to be satisfied. ) Near p = 0, cj){p) - cj){Q) -f Cp^ + Filnp, with 
C = |S|if'(0(O),6')/4. The requirement Fi = and D2 = Q determines 0(0) and Di. E is fixed by the 
total flux 2tt j pdpe{(p{p))Ez = g. Now we look for the solution which increases monotonically from 0(0) 
to 0(00) = 0c. This requires C > 0, which is shown to be consistent afterwards. In order for such a 
solution to exist, -ftr(0, 9) should have a bump in the region < < 0c. To see this, we note that 0"(p) 
changes sign at some pd, where (j)"{pd) = 0. From and <p'[pd) > 0, K'{(j){pd),9) has to be positive, 
which connects smoothly to the behavior near cj) = (j)c, where K' < 0. So K' = at some = 0^, the 
bump position, in the range (pd < 4>s < 4>c- If we define ps by 4>{ps) = 0s, it is the measure of the radius 
of the tube. Below, > 1 is proved for any solution, so we need the bump structure for 9 > 1. This 
condition is met by the existence of the bump for 6* = 1, as is clear from K{{), 9) ~ —9. (see the Figure. 
) Thus we can say that the stability of the condensed vacuum assures the flux tube of infinite length. 

To show 9 > \, 'A sum rule is derived by multiplying d(j)/dp on both sides of ([9]) and integrating from 
p = to 00. Using boundary conditions on both ends, we get 

1 f d^\^ ^ / ^^^^ _ ^^^^^^^ ^ 



P\ dp J 



Since K{(t)c,d) ~ —1, K{(t){0),9) < —1 holds, implying 6 > 1; the gain of the electric energy inside the 
tube is larger than the loss in breaking the condensation energy. If the left-hand side of ([9]) is regarded as 
the surface anergy, the difference of these two energies is sustained by the surface energy. If we neglect 
completely the surface energy, the solution satisfies K'{(f),6) — 0, thus (p = (f)c, or = 0. Setting = 
inside the tube and cj) — (pc outside, and denoting by S the cross section of the tube, the energy per 
unit length S{\B\ + E"^ /2) is minimized under the constraint SE = g. Then S = .g/^2|B|, E = \/2\B\ 
is obtained, leading to 9 = I. The same result of course can be derived by the sum rule since it gives 
K{(p{0), 6) = —1 and </)(0) = if the surface energy is neglected. These hold for all a > 1. If the surface 
energy is included and becomes large, 9 and hence E increases {E < -Emax. when a — 1) and both the 
bump position and (/'(O) approach (j)c. Thus a solution to ([9|) is assured again for any a > 1. 

Near the source; 

In the neighbourhood of the point-like colored souce of a quark, the condition ([4]) near cf) — becomes 
essential. The equation to be solved is = V'{(j)) — e'(0)i?"^/2. We assume (j) approaches near 

the source and search for the spherical solution with the electric field at the distant r from the quark 
given by E"^ = g\^/2r^. Neglecting V'{(l)) and by e{(j)) ~ 1 + acj)'^ (a < 0) near cj) — 0, we get 4>{r) = 
(l/r) exp(— \/A/r), A = —ag^^^iy^ /2)^ ~ > 0. The solution is consistent with the starting 

assumption that ~ near the source. The effect of the condensation near the quark is smaller than 
any power of r. 

Strong coupling for the canonical field; 

Our full quantum theory is defined by regarding ([l]) as the action which is functionally integrated 
over (p and v4°(a;). Let us neglect the fluctuation of cj) hi'st. In the condensed vacuum, since e((/)c) — 0, 
fluctuate indefinitely but it has no observable effect by the same fact e = 0. In the homogeneous 
case, Ce^A can be written by the canonical field ^^(a;) — y/e{(t))A'J;^{x) (unity for the coefficient of the 
kinetic term). Then we get the same Lagrangian as QCD with g replaced by g^ — gj \/ e(0))- Thus 
QCD becomes a strong coupling theory near the vacuum </) — 0c, where g^ — s- oo. Then, the dominant 
interaction term becomes 4-point vertex and by the stationary phase machanism, all ^^(x)'s are freezed 
to zero locally; gluons are confined in the operator level. The mechnisms of confinement of gluons are 
different for and but they are physically equivalent. When the quark field q{x) are introduced, 
we add the usual term q(a;)(7^(i9p + gA^J_(x)) -\- m)q{x) to CcS with = X]a=i ^Ij,^"" /'^- Note that the 
gluonic condensation is unrelated to the quark sector. When e = 0, the integration by A'^{x) leads to 
the vanishing of the operator local color current of quarks; j^{x) — q{x)jn{X°' /2)q{x) — for any a and 
/i. To study excitations, we note that fluctuations of A^ and come about only if (f> fluctuates from the 
value (f>c, which is governed by the glue ball mass. We expect (f> ^ (f>c in a localized region with size \/mg, 
where A^'s and have non-zero values and the ordinary perturbative picture works. Above picture of 
strong coupling expansion looks very similar to the same expansion in lattice QCD. 

Discussions; 



We try to interpret our results by the dual terminology. The permeability /i = 1/e is infnite in 
the vacuum and for the above tube-like solution, the magnatic supercurrent defined by rotI?°/e = 
(l/e)Ve X £)" is circulating around the surface of the tube which sustains the electric flux. Let us 
introduce the dual potential as a Lagrange multiplier of the Bianchi identity [TH [T3] by taking the the 
axial gauge n^A"^ = with some constant vector n^. By a simple manipulation, we get the relation that 
rotC° is equal to plus the string term which depends on Thus the dual potential C, with the 
string part subtracted, constitutes the tube decaying to zero for p — oo. Although we have not explicitly 
the mass term rn^C^ which breaks the magnetic gauge invariance, the model discussed here behaves as 
dual to the superconductor in the above sense. 

Our starting formula including full quantum effcts is ([8|), which can be used in principle to determine 
consistently phenomenological parameters introduced above, or to study the shielding charge due to 
gluons, or by including quark degrees, the quark pair production which breaks the tube, e.t.c. 

References 

[1] G. B. Bah, C. Schhchter and K. Schilling, Phys. Rev. D51, 5165 (1995); R. W. Haymaker, V. Sihgh, 
Y. Peng and J. Wosiek, Phys. Rev. D53, 389 (1996); A. M. Green, C. Michael and P. S. Spencer, 
Phys. Rev. D55, 1216 (1997); Lattice 2004, Proc. XXIInd Int.Conf.on Lattice Field Theory, Batavia, 
IL, USA, June 2004; M. A. Stephanov, arXiv : hep-lat /0701002| v^l. 

[2] G. 't Hooft, Nucl. Phys. B190 [FS3], 455 (1981). 

[3] H. Toki and H. Sugawara, Prog. Part. Nucl. Phts. 45, S397 (2000); A. Di Giacomo, Nucl. Phys. A 
702, 73c (2002); arXiv:hep-lat/0710.1164vl; J. Greensite, Prog. Part. Nucl. Phys. 51, 1 (2003) G. 
Ripka, Dual Superconductor Models of Color Confinement^ Lect. Note Phys. 639 ( Springer- Verlag 
Berlin Heidelberg 2004), DOI 10.1007/b94800; V. G. Bornyakov, M. I. Pohkarpov, G. Schierholz, 
T. Suzuki and S. N. Syritsyn, Nucl. Phys. Proc. Suppl. 153 25 (2006); K-I. Kondo, A. Shibata, T. 
Shinohara, T. Murakami, S. Kato and S. Ito, larXiv:0803.2451K ^l [hep-lat]. 

[4] V. N. Gribov, Nucl. Phys. B139, 1 (1978); D. Zwanziger, Nucl. Phys. B518, 237 (1998); Y. Nakagawa, 
A. Nakamura, T. Saro and H. Toki, Phys. Rev. D75, 104508 (2007). 

[5] J. Kogut and L. Susskind, Phys. Rev. D9, 3501 (1974). 

[6] G. 't Hooft, Quarks and gauge fields, in Proc. CoUoq. on Lagrangian field theories 
(C.N.R.S.,Marseille,1974) 

[7] M. A. Schifman, A. I. Vainstein and V. I. Zakharov, Nucl. Phys. B147, 385, 448, 519 (1979); R. 
Fukuda and Y. Kazama, Phys. Rev. Lett. 45, 1142 (1980). 

[8] A. M. Polyakov, JETP Lett, 20, 194 (1974) [Pisma Zh. Eksp. Teor. Fiz. 20, 430 (1974)] ; G. 't Hooft, 
Nucl. Phys. B 79, 276 (1974). 

[9] R. Fukuda and T. Kugo, Prog. Thor. Phys. 60, 565 (1978); R. Fukuda, Phys. Lett. 73B, 305 (1978); 
60, 1152 (1978); Phys. Rev. D21, 485 (1980). 



[10] T. Kugo and I. Ojima, Prog. Theor. Phys. Suppl. 66, 1 (1979). 
[11] M. Gell-mann and F. Low, Phys. Rev. 84, 350 (1951). 
[12] M. B. Halpern, Phys. Rev. D19, 517 (1979). 

[13] R. Fukuda, Prog. Theor. Phys. 67, 648 (1982); 67, 655 (1982); 68, 602 (1982). 



